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7": ABSTRACT 

u 

is^ ; Recent studies of type la Supernovae with redshifts up to about z < 1 reveal evidence for 

H ' a cosmic acceleration in the expansion of the Universe. The most straightforward explanation 
to account for this acceleration is a cosmological constant dominating the recent history of our 
Universe; however, a more interesting suggestion is to consider an evolving vacuum energy. 
Several proposals have been put forward along these lines, most of them in the context of 
General Relativity. In this work we analyse the conditions under which the dynamics of a 
self-interacting Brans-Dicke field can account for this accelerated expansion of the Universe. 
We show that accelerated expanding solutions can be achieved with a quadratic self-coupling 
of the Brans-Dicke field and a negative coupling constant uj. 
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1 Introduction 

Type la Supernovae (SNe la) allow us, when used as standard candles, to reliably determine 
the cosmological parameters. Indeed, low redshift supernovae can be used to obtain the Hubble 
constant, Hq, while supernovae at greater distances allow probing the deceleration parameter, 
go- Recent analyses of the magnitude-redshift relation of about 50 SNe la with redshifts greater 
than z > 0.35 strongly suggest that we are living in an accelerating, low-matter density Universe 
[1 - 4] . The consistency relationship between these cosmological parameters and the luminosity 
distance, d^, of a SNe is given, for z < 1, hj d^ ^ Hq-^Iz + (1 — qo)z'^/2]. The results obtained 
by different groups indicate that [1 - 4] 

- 1 < go < . (1) 

These values seem to favour, in the case of a flat Universe, a rather low contribution to the 
energy density from non-relativistic matter, say Qm ~ 0.3, but on the other hand a dominant 
positive cosmological constant, Q\ ~ 0.7, as can be inferred from data and consistently from 
the well known Friedmann model equation 

qo = ^{S-f + 1)Qm - Qa , (2) 

where fiM{v) denotes the energy density of matter (vacuum) in terms of the critical density, 7 
stands for the constant in the matter equation of state, p = 7p, and observationally — 1 < 7 < 
(7 = for non-relativistic matter, 7 = 1/3 for relativistic matter and 7 = —1 for a scalar field 
fluid dominated by its potential). Even though a dominant cosmological constant corresponding 
to Q/y ~ 0.7 is consistent with all observational facts (see for a list of constraints and for a 
thorough discussion) it implies a quite unnatural fine tuning of parameters if its origin lies in any 
known particle physics scale (see Refs. [^ for reviews and Refs. |^, ^ for possible connections 
with fundamental symmetries such as Lorentz invariance and string theory S-duality). 

From Eq. (^ we see that a negative go can also arise from a negative pressure component or 



"dark energy" . This alternative calls for a scalar field endowed with a potential which can give 
rise to a dynamical vacuum energy, the so-called "quintessence" . Suggestions along these lines 
were proposed long ago, although yielding, in that case, a vanishing deceleration parameter 
[p!0| ]. A number of quintessence models have been put forward, most of them invoking a scalar 
field with a very shallow potential, which until recently was overdamped by the expansion 
of the Universe, so that its energy density was smaller than the radiation energy density at 
early times 0, 12, 13|. For instance, scalar fields with an exponential type potential can. 



under conditions, render a negative go iHl [Hi- More involved possibilities include the string 



theory dilaton together with gaugino condensation [|T6|, an axion with an almost massless 
quark |T^, the effect of i^-particle recoil |]TB[, supersymmetric inflationary models |TP|, the 



multidimensional Einstein- Yang-Mills system pQ], etc. However interesting, most of these 



suggestions necessarily involve a quite severe fine tuning of parameters |21]. This fact calls for 
constructions that allow for a negative deceleration parameter using sources of negative pressure 
that do not require a potential arising from known particle physics scales. Other difficulties 
associated with quintessence is that the couplings of the scalar field to matter can lead to 
observable long-range forces and time variation of fundamental constants of nature [E^, I2H] . 



In this work we study the possibility of obtaining the required negative pressure effects by 
considering a scalar field coupled non-minimally with gravity. It is well known that the non- 
minimal coupling of fields to gravity have non-trivial implications on issues such as geodesic 
completeness, stability of the ground state (see p4| and references therein) and the variation 



of constants of nature. In here, we shall consider the self-interacting Brans-Dicke theory as a 
prototype of a non-minimally coupled theory and show that accelerated expanding solutions 
can be obtained and the abovementioned difficulties partially evaded. Actually, accelerated 
expanding solutions at present time, the so-called scaling attractors, have already been studied 
in some theories of gravity with non-minimal coupling to a scalar field. Indeed, the non- 
minimal coupling \e(tP'R has been considered in Ref. |^ with the conclusion that the scalar 



field behaves as a barotropic fiuid, then leading to scaling attractors, only when the constant 



e << 1 for exponential and power-law potentials. In Ref. [^] the search of scaling attractors 



has been extended to the case of the non-minimal couphng of the form [1 + e/(0)]i? and 
V{(f>) = yl/(0)*^, where f{(f>) is a power-law or exponential function of the scalar field and 
A and M are constants, and shown that, in the limit where the kinetic term of the scalar 
field is dynamically unimportant, the constraint on the time variation of the gravitational 
coupling is quite severe and limits the fraction of energy in the scalar field to be at most 
4% of the total energy density. We shall see in the next sections that our results, although 
obtained from somewhat different starting assumptions, are consistent with the conclusions 
of Refs. ||2^, |2^, namely that although accelerated expanding solutions can be obtained in 
non-minimally coupled gravity theories the constraint on the variability of the gravitational 
coupling is quite strong and that it implies a Universe that is considerably older than Hq^. 

Our starting point is the Brans-Dicke (BD) theory [^ as this is a viable scalar-tensor 
alternative to General Relativity. In this theory spacetime is described by the interplay of the 
metric tensor with a scalar field, 0, so that the strength of the gravitational coupling to matter 
is given by 0~^. The theory is consistent with observational tests which explains the renewed 
interest in its application both in astrophysics and cosmology. Moreover, the BD theory is also 
the model arising from string theory at low-energies in the so-called string or Jordan frame and 
from the dimensional reduction of Kaluza-Klein theories. BD theories and its extensions are 



known to have relevant cosmological implications [^ ^, [3^ . 



In its simplest version, the salient features of BD theory depend uniquely on the strength 
of the dimensionless "Dicke coupling constant", u, that couples the scalar field universally to 
matter. However, the cosmological setting arising from the simplest BD scenario is, as far 
as power-law solutions are concerned (cf. equations, below), inconsistent with the cosmic 
acceleration, unless \u!\ is of order unity. Since the coupling u is observationally constrained by 



solar-system experiments to be |ci;| > 500 [^, we are led to consider a version of the BD theory 
where the BD scalar has a potential. For the latter we choose a quadratic self-coupling in the 
Jordan frame as this implies in a minimal change in the field equations. Thus, our scenario can 
be regarded as the next to minimal Brans-Dicke model. Furthermore, we show that a negative 
uj is required to account for the accelerated expansion of the Universe. 



In what follows we consider an homogeneous and isotropic Friedmann- Robert son- Walker 
Universe where matter (not including the BD scalar field) is described as a perfect fiuid. The 
relevant thermodynamical fluid variables are the energy density and the isotropic pressure that 
are related by the equation of state specified above. 

2 Field Equations and Solutions 

The self-interacting BD field theory of our interest is described by the following action: 



S 



16tt 



j4 / 

a xy/—g 
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+ V{(P) 



+ s. 



Matter 



(3) 



where we choose V{(f)) = Vo0^. As already mentioned in the Introduction a BD theory with a 
vanishing potential is incompatible with a negative go; at the least for power-law solutions for 
the scale-factor and scalar field specified below, unless |u;| = 0(1) which is inconsistent with 
solar system experiments. However, as we shall see in a while, the potential we have introduced 
for the scalar field changes the equations of motion in a minimal way. It is worth stressing 
that introducing a potential does not involve, in our proposal, any fine tuning of parameters 
as our scale, Vq, will be fixed by the value of Newton's constant and the age of the Universe. 
Furthermore, as far as giving origin to a time variation of the fundamental constants, our 
proposal implies indeed in a sizable time variation of the gravitational coupling, even though 
this is compatible with available data. 

From the Lagrangian density (H) we obtain the field equations: 
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(5) 



Hence for our choice of the potential, V{(j)), the last term in Eq. (^ vanishes and V{(f)) affects 
directly only the dynamics of the scale-factor, a{t). 

These field equations in a Friedmann-Robertson- Walker geometry 



read 



ds^ = -dt^ + a{tf 



dr"^ 
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Clearly, these equations must be considered together with the energy conservation equation 



p + 3-{p + 3p)=0 
a 



(10) 



arising from the covariant conservation of the energy-momentum tensor of matter that has been 
assumed to behave as a perfect fluid. 



Finally, the relationship between the gravitational coupling G{t) and the scalar field is given 



by 



^^^^ - [^rd W) 



(11) 



In order to obtain solutions for the above equations we consider the following power-law 
form for both scale-factor and scalar field: 



a{t) = ao{ — 

.to 



t\^ 



0(t)^0o(-) 



(12) 



where the zero indices stand for the present time. 



Substituting these solutions into Eq. (||) leads, when neglecting the mattter pressure con- 
stribution (7 ~ 0), to the relationships 



3a = 2-13 ; /3 



a 



3 + 2uj 



(13) 



where 



a = 



Snpotl 



(14) 



From Eqs. (|ri 



- ([T3|) and the Hubble parameter at present, Hq = ^{to), we get the deceleration 
parameter, go = ~^(^o)) and the time variation of the gravitational coupling: 



?o 



l + P 
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Moreover, adjusting the gravitational coupling to its present value, Gn, i.e. Newton's constant, 
we get, after inserting Eq. (|Tl|) into Eq. (|l^ and combining with Eq. (|13|), for the age of the 
Universe 



t 



2P{iu + 2) 



3QmHq 

Naturally, for consistency Vq must satisfy Eq. (|^). Hence, for a flat Universe {k = 0) we find 



17) 



Vo = 3Gn^mHI{2uo + 3) 



2(3{2u + 3) - 6a(a + /?) + uf3'^ 
4/?(cu + 2)2 



3 Conditions for negative qq 



In order to have go < it is equivalent that the right hand side of Eq. (§) is positive. It is 
useful then to consider it as a function 



/--(3^i'--)^i-Kff-^:-y^ <-) 



and establish the conditions for having / > 0. 



Power- law solutions Eq. (|T^) satisfy the field equations only if /3 = —2. Then, solving for lu 
yields 

.<-y^-i. (20) 

Since the product Vq^o^o i^ completely determined by Eqs. (Illl), ([T7|) and ([T^) it is easy to see 
that this condition is satisfied for any value of uj. Therefore, we can conclude that the cosmic 
accelerated expansion can be driven by a self-interacting Brans-Dicke field as we have specified. 

In summary we have, after establishing that (3 = —2 



ait)=ao[^y ; 0(t) = 0o (^) , (21) 



^2_ Aju + 2) 
3QmHq 



^0 - oo E7-2 ' y'^'^) 



and 



l)r^"°' '"' 



Qo = —. (24) 



We see that the obtained value for go is consistent with the observations, Eq. ([I|), and so is 
the resulting time variation of the gravitational coupling [p^, even though many searches are 



consistent with no variation at all |32]. We stress that, in our proposal, the gravitational cou- 
pling is a growing function of time. Morever, notice that, although / is positive independently 
of the sign of cu, as (3 is negative then in order to have a meaningful age of the Universe it 
implies that the Dicke coupling constant must be negative. Interestingly, negative values for 
uj are found in the BD effective low-energy models arising from Kaluza-Klein and superstring 



theories |0. We shall see in the next section that negative values for uj are also required for 



obtaining growing modes for the energy density perturbations in a Universe expanding in an 
accelerated way. 

4 Perturbed Field Equations and Asymptotic Behaviour 

Having found accelerated expanding solutions in the context of the BD theory we should ask 
whether the issue of structure formation is modified by the dynamics of the BD scalar field. 
That is, we are bound to consider the evolution of the energy density perturbations in the 
context of our BD scenario as these are, of course, associated with the formation of structure 
in the Universe. Hence in what follows we shall develop a formalism that quickly allows us to 
find the asymptotic behaviour of the relevant variables for accelerating expanding solutions at 
present time. This treatment is similar, in its results, to the more encompassing analysis of 



cosmological perturbations for generalized theories of gravity ||3^ , at least in what concerns the 
asjTiiptotic behaviour of solutions. 

In order to get the matter energy density perturbations we consider the temporal compo- 
nents of Eqs. (H), (P) and ([T0|) after perturbation. Then to obtain the relevant perturbed 
equations we write the metric tensor as 

9tiv -^ gfiu + ^Qi^y , (25) 



where 5g^y = h^^ , and work in the gauge /iq^ = 0. Our conclusions will be independent of this 



gauge choice (see ||33| for a discussion). For this perturbed metric [^ 



5-Rnn = ^T 



hkk — "^-hkk + 2 
a 



Ikk 



The perturbation of the energy-momentum tensor is given by 



(26) 



5T^" = 5p 



(27) 



and the corresponding trace reads 



ST = Sp — 3Sp . 



For the perturbation of the d'Alembertian of the BD field we have 



1_: 






a a^ 



The relevant perturbations are parametrized in the following way: 



(2^ 



(29) 



hkk — ah 



(30) 



S(j) = Xcj) , A << 1 



(31) 



Sp = Ap , A « 1 , (32) 

where h{t), X{t) and A(t) are the perturbed gravitational field, scalar field and matter energy 
density, respectively. 
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Since structure is formed when pressure no longer prevents gravitational collapse we set 
7 = 0. The perturbed equations are then 



^'-D'"-^|:^'--'-^^-^'-4^-T^^' '-) 



^+H^4)^-(^«^|)-^|'-^-l3l^-' <-) 



A-^h + 6U', = , (35) 

where Un is the comoving fluid velocity. 



Inserting solutions (0), (|1^) with /5 = — 2 into this set of coupled differential equations we get 



lU^/. = /?(2 + c.)(^H-A + 2(l + .)fA-f0o(^)'A 






A-i/i + <5f/t = . (38) 

In order to continue we suppose that the perturbations behave as plane waves: 

A(f, t) = X{t)exp{-ik ■ r) , (39) 

where k is the wave number of the perturbation, and we set 5U^ to vanish, which is allowed by 
an infinitesimal gauge transformation. 

Combining then Eq. (p6D with (|38|) we get 
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A.(1)(A-.).(^)(A-A)., , 






A + -^^^T^A + ^(^) A = . (40) 



Since we are interested in the asymptotic regime we keep only terms up to t ^ and neglect the 
last term of the previous equation. 

Finally, from the value of the product, Vq^o^o — ~(4^ + ^) previously obtained, we get the 
differential equation 



where Ci = y. 



A + ^(A-A) + ^(A-A) + ^A + ^A = 



(41) 



Ciioj) =Q + 2u 



(42) 



C^{uo) = A{1 + uj) + Ci 



(43) 



and 



CAu) = 



lO + Qio 



(44) 



Aiming to solve Eq. (^if) we look for solutions of the following form: 



A-A = /(t) ; f{t)=^t' ; A{t)=xt' 



(45) 



where ^ and x are constants. 



Inserting these solutions into Eq. ( |1T| ) we obtain the algebraic equation 



XO' + [x{Cs - 1) + e(Ci - Cs)] 9 + Cix -O+^C2 = 



(46) 
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and also that 6 = 6. 



Solving for 6 we find: 



. X(l - Cs) + aCs - Ci) ± V[X(C3 - 1) + aCi - Cs)r - Ax [C,ix - + ^C,] 

^± = 7^ ■ [Ai) 



2x 

Thus, the existence of growing modes for the energy density perturbations corresponds now 
to whether there are positive 6 solutions. In the case of our interest, namely u negative and 
\uj\ >> 1, it is easy to see that 



2\uj\ 



-«*, 






2\lj\ 



2^ 

Xj 



(4J 



meaning that 6*+ corresponds to a growing mode whenever ^ < X- Moreover, we see that 
it is only for u negative that real 6 solutions always exist. We can then conclude that the 
cosmological setting of our accelerated expanding solutions does not upset known structure 
formation scenarios. 

5 Discussion and Conclusions 



In this work we have shown that a minimally extended BD theory with a quadratic self-coupling 
in the Jordan frame and a negative uj can account for the accelerated expansion of the Universe 
yielding go = ~\- The resulting variation on time of the gravitational coupling is given by 
(^ j = I iJo) being still compatible with data. In this respect our results are similar to the 
conclusions of Refs. [^ |26| as likewise we find that the existence of scaling attractors in 
gravity theories with non-minimal coupling to a scalar fields is severly constrained by the time 
variability of the gravitational coupling. Furthermore, we have also shown that the model allows 
growing modes for the energy density perturbations of matter implying that the dynamics of 
the BD field does not upset known structure formation scenarios . 
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Before closing, it is worth pointing out some distinct features of our model. The first 
one being that as accelerated expanding solutions require a negative u;, we then expect the 
parametrized-post-Newtonian (PPN) parameter •Jppn = ^rf; to be, for large |ti;|, fairly close 
but greater than 1. This is a clear observational signature of our proposal, which even though 



being consistent with current data [p5|] 



\-fPPN-l\ < 2 X 10~^ , (49) 

can be, at least in principle, further improved by dedicated solar system experiments. 

A second consequence of our proposal is that it implies the Universe is considerably older 
than Hq^. This is incompatible with observation if the age of the Universe is identified with 
age of the oldest stars, the globular clusters. Contradiction with models of chemical evolution 
of galaxies may also exist. Even though these constraints are extracted from data assuming 
a Universe with constant gravitational coupling they point out that the minimal extention 
of the BD model we have considered is not quite consistent with data. In this respect, it is 
interesting that the effect of variation on time of the gravitational coupling may reveal itself in 
the evolution of astrophysical objects. Indeed, as recently discussed [Q, stringent bounds on uj 
can be set, namely that cu > 5000 and that f ^j = O(10~^^), based on the luminosity function 
of white dwarfs assuming that G < and 12.5 Gyr as the age of the Universe. We suspect that 
this age of the Universe problem is a common feature of all scalar-tensor gravity models that 



have scaling attractors and do not contain, as for instance in p6|, the Einstein-Hilbert term in 
the action. 

Finally, we could conclude remarking that an interesting theoretical challenge would be de- 
vicing inflationary models which would, at late times, behave like scalar-tensor gravity theories 
in what concerns the existence of scaling attractors. 
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